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O ! 

, In the framework of the fci-factorization approach, the production of T mesons at 

^ ' the Fermilab Tevatron and CERN LHC is considered, and the predictions on the spin 

Hi ! 

■ alignment parameter a are presented. We argue that measuring the polarization of 

quarkonium states can serve as a crucial test discriminating two competing theoretical 

^ ' approaches to parton dynamics in QCD. 

ft, 

I. INTRODUCTION 

> ' 

Nowadays, the production of heavy quarkonium states at high energies is under intense 

(N ■ . n n 

theoretical and experimental study [2|J. The production mechanism involves the physics 
of both short and long distances, and so, appeals to both perturbative and nonperturbative 
methods of QCD. The creation of a heavy quark pair QQ proceeds via the photon-gluon or 
gluon-gluon fusion (respectively, in ep and pp collisions) referring to small distances of the 
order of l/(2mg), while the formation of the colorless final state refers to longer distances of 
the order of l/[niQa s (mQ)}. These distances are longer than the distances typical for hard 
interaction but are yet shorter than the ones responsible for hadronization (or confinement). 
Consequently, the production of heavy quarkonium states is under control of perturbative 
QCD but, on the other hand, is succeeded by nonperturbative emission of soft gluons. This 
feature gives rise to two competing theoretical approaches known in the literature as the 
color-singlet [j], Q| and color-octet p] models. According to the color-singlet approach, the 
formation of a colorless final state takes place already at the level of the hard partonic 
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subprocess (which includes the emission of hard gluons when necessary). In the color-octet 
model, also known as nonrelativistic QCD (NRQCD), the formation of a meson starts from a 
color-octet QQ pair and proceeds via the emission of soft nonperturbative gluons. The former 
model has a well defined applicability range and has already demonstrated its predictive 
power in describing the J/tf> production at HERA, both in the collinear p3|_ and the k t - 
factorization [3] approaches. As it was shown in the analysis of recent ZEUS Q| data, there 
is no need in the color-octet contribution, neither in the collinear nor in the fc^-factorization 
approach. 

Originally, the color-octet model was introduced to overcome the discrepancy between the 
large J/ if> production cross section measured in pp interactions at the Tevatron p, LLC), and 
the results of theoretical calculations based on the standard perturbative QCD technique. 
The problem was apparently solved by attributing the discrepancy to the hypothetical 
contributions from the intermediate color-octet states, which must obey certain hierarchy 
in powers of the relative velicity of the quarks in a bound system. However, the numerical 
estimates of these contributions extracted from the analysis of Tevatron data are at odds 
with the HERA data, especially as far as the inelasticity parameter z = E^/E 1 is concerned 



12|. In the fc t -factorization approach, the values of the color-octet contributions obtained as 



fits of the Tevatron data appear to be substantia 
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ly smaller than the ones in the collinear 



scheme, or even can be neglected at all Jjl, 

In the present note we want to stress once again that measuring the polarizaton of 
quarkonium states produced at high energies may serve as an important and crucial test 
discriminating the different theoretical concepts. The first attempts to solve the quarkonium 
polarization problem within the A; r factorization approach were made in the pioneering work 
17] (see also for ep collisions and in Refs. [lj, 111]] for pp collisions. It was emphasised 



that the off-shellness of the initial gluons, the intrinsic feature of the fc r factorization 
approach, has an immediate consequence (by analogy with longitudinal photons) in the 
longitudinal polarization of the final state J/ip mesons. The theoretical predictions [l7|] have 
stimulated experimental investigation of J /if} spin alignment at the collider HERA. The 
first results obtained by the collaborations HI and ZEUS have been described in Ref.{3|. 
These results have qualitatively confirmed the predictions on the dominance of longitudinal 
polarization. 

The preliminary results on the J /if) polarization at the Tevatron obtained by the 



collaborations E537 [20] and CDF [2lJ also point to logitudinal polarization with the average 
value of spin alignment parameter a ~ —0.2 over the whole range of J/ip transverse 
momentum px- The collaboration DO is currently analysing the data on the spin alignment 
of T mesons. 

In the NRQCD approach, the problem of quarkonium polarization remains unsolved 
l|,l2j. The gluon fragmentation mechanism leads to strong transverse polarization. Including 
the next-to-leading QCD corrections makes the transverse polarization even stronger. The 
only way out is seen in increasing the fraction of unpolarized mesons by attributing large 
contributions to the certain color-octet channels which, however, violates the expected 
NRQCD hierarchy. The role of the color-octet contributions taken into account in the 
analysis of recent ZEUS data is obscure and does not lead to a conclusive description of 
J/ip polarization parameters (see {2I). 

II. NUMERICAL RESULTS 

The goal of this paper is to derive theoretical predictions on the polarization of T mesons 
produced at the Fermilab Tevatron and CERN LHC. In the fc r factorization approach, the 
cross section of a physical process is calculated as a convolution of the partonic cross 
section a and the unintegrated parton distribustion T g {x, k^, /1 2 ), which depend on both 
the longitudinal momentum fraction x and transverse momentum kr' 

CTpp = J J~ g {X\ , k^rp } jj, ) J- g (x2 1 ^2T> ) ® 99 (^'l ; ^2 1 ^2Ti " •) ^X\ d,X2 dk-^j> rfA^j'. (1) 



In accord with the fc t -factorization prescriptions [22|, |23j, |2J, [25|], the off-shell gluon spin 
density matrix is taken in the form 



^ = ^y pX y\k T \ 2 = e T k v T /\k T \\ (2) 

In all other respects, our calculations follow the standard Feynman rules. 

In order to estimate the degree of theoretical uncertainty connected with the choice of 



unintegrated gluon density, we use two different parametrizations, which are known to s 



now 



the largest difference with each other, namely, the ones proposed in Refs. 



In the first case 



22 



25| and 



26|]. 



22|, the unintegrated gluon density is derived from the ordinary (collinear) 



density G(x, /z 2 ) by differentiating it with respect to /i 2 and setting /i 2 = k\. Here we use the 
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LO GRV set [27| as the input colinear density. In the following, this will be referred to as 



dGRV parametrisation. The other unintegrated gluon density 



of leading order BFKL equation 



26j] is obtained solution 



25| in the double-logarithm approximation. Technically, it 



is calculated as the convolution of the ordinary gluon density with some universal weight 
factor. In the following, this will be referred to as JB parametrisation. 

The production of T mesons in pp collisions can proceed via either direct gluon-gluon 
fusion or the production of P-wave states Xb followed by their radiative decays Xb~ >T+7. 
The direct mechanism corresponds to the partonic subprocess 



g + g -> ^ + g 



(3) 



which includes the emission of an additional hard gluon in the final state. The production 
of P-wave mesons is given by 

9 + 9 ^Xb, (4) 

and there is no emittion of any additional gluons. As we have already mentioned above, we 
see no need in taking the color-octet contributions into consideration. 

The other essential parameters were taken as follows: the 6-quark mass m b = m T /2 = 4.75 
GeV; the T meson wave function |\l/x(0)| 2 = 0.4 GeV 3 (known from the leptonic decay 



width T 



i+i- 



28j); the wave function of P-wave states (0)| 2 = 0.12 GeV 5 (taken from the 



potential model [29|]); the radiative decay branchings Br(xb,j^^l) = 0.06, 0.35, 0.22 for 
(J = 0, 1, 2) [281] ; the renormalization and factorization scale /i 2 = m\ + p\. 



The polarization state of a vector meson is characterized by the spin alignment parameter 
a which is defined as a function of any kinematic variable as 



a (V) = (da/dV - 3da L /dV)/(da/dV + da L /dV), 



(5) 



where a is the reaction cross section and is the part of cross section corresponding to 
mesons with longitudinal polarization (zero helicity state). The limiting values a = 1 and 
a = — 1 refer to the totally transverse and totally longitudinal polarizations. We will be 
interested in the behavior of a as a function of the T transverse momentum: V = \pr\- The 
experimental definition of a is based on measuring the angular distributions of the decay 
leptons 

dT(T^fi + fi-)/dcose ~ 1 +ctcos 2 #, (6) 
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where 6 is the polar angle of the final state muon measured in the decaying meson rest 
frame. 

The results of our calculations for the kinematic conditions of the Tevatron and LHC 
are displayed in Figs. 1 and 2. In both cases, the integration limits over rapidity were 
adjusted to the experimental acceptances of CDF (\y?\ < 0.6) at the Tevatron and ATLAS 
(|yx| < 2.5) at the LHC. The upper panels show the predicted transverse momentum 
distributions. Separately shown are the contributions from the direct (dashed lines) and 
P-wave decay (dotted lines) mechanisms. Note that, in spite of the suppression due to 
smaller values of P-wave finctions compared to S'-wave finctions, the dominant contribution 
comes from the subprocess ([4j) rather than from j3]). The reason can be seen in the much 
smaller values of the final state invariant masses, m x « m^ g . Our conclusion on the 
relative size of two contributions is compatible with the preliminary estimates obtained by 



the collaboration CDF 



The pt shape of the individual contributions yet has not been 



meaured experimentally. In the fc t -factorization approach, this shape is determined by the 
unintegrated gluon density. The average px is a bit lower in the process (J3|), because the 
total transverse momentum (equal to that of the initial gluons) is shared between the two 
final state particles; at the same time, the contribution from the matrix element is nearly 
unimportant. 

It is worth noting that the production of Xb mesons can hardly be described in a consistent 
way within the collinear factorization scheme. The leading order contribution coming from 
the subprocess ([4]) shows unphysical 5-like px spectrum. The usual lame excuses that the 
particles produced at zero px disappear in the beam pipe and remain invisible do not work, 
because the decay products do have nonzero px and, undoubtedly, can be detected. At the 
same time, the introduction of next-to-leading contributions (i.e., the processes with extra 
gluons in the final state) causes the problem of infrared divergences. 

The central panels in Figs. 1 and 2 show the behavior of the spin alignment parameter 
a for T mesons produced in the direct subprocess j3]). The increase in the fraction of 
longitudinally polarised mesons is promptly connected with the increasing virtuality (and, 
consequently, the strenghtening longitudinal polarization) of the initial gluons. 

As far as the decays of P-wave states are concerned, nothing is known on the polarisation 
properties of these decays. If we assume that the quark spin is conserved in radiative 
transitions, and the emission of a photon only changes the quark orbital momentum (as it 
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is known to be true in the electric dipole transitions in atomic physics, AS = 0, AL = ±1), 
then the predictions on a appear to be similar to those made for the direct channel (see lower 
panels in Figs. 1 and 2, dotted curves). If, on the contrary, we assume that the the transition 
Xb^T + 7 leads to complete depolarization, then we arrive at a more moderate behavior 
of the parameter a (dash-dotted curves in Figs. 1 and 2). The overall polarization remains 
slightly longitudinal (a ~ —0.2) in the whole range of pr due to the 'direct' contribution. A 
comparison between the data on J / ip and ip' polarization at the Tevatron [21fl seems to give 
support to the depolarization hypothesis. The difference between the J/ip and ip' polarization 
data can be naturally explained by the presence of the depolarizing contribution in the case 
of J /ip and the absense of this contribution in the case of ip' . 

A state with purely direct production mechanism in the bottomonium family is the T(3S f ) 
meson. The calculations presented here are also valid for this state, except the lower total 
cross section (by an approximate factor of 1/3) because of the correspondingly lower value 
of the wave function |\I/t(3S)(0)| 2 = 0.13 GeV 3 . At the same time, the predictions on the spin 
alignment parameter a remain intact (central panels in Figs. 1 and 2). 



III. CONCLUSIONS 



We have considered the production of T mesons in high energy pp collisions in the kt- 
factorization approach and derived predictions on the spin alignment parameter a(pr). We 
point out that the predicted value of a(pr) is typically negative in the whole range of px 
and shows variations from a ~ (—0.2) to a ~ (—0.7) depending on the hypothesis assumed 
for the decays Xb~ ►T(lS')+7. At the LHC energies, the theoretical predictions possess less 
sensitivity to the choice of unintegrated gluon distributions. The purest probe is provided 
by the polarization of T(35*) mesons. In that case, the polarization is the strongest and the 
predictions are free from uncertainties coming from the unknown properties of Xb decays. 

We do not discuss the behavior of the parameter a(px) at asymptotically large transverse 
momenta where the applicability of the /c t -factorization approach is questionable. 
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Phc. 1: Predictions on the production of T mesons at the Tevatron. Thick lines, JB 
parametrization; thin lines, dGRV parametrization. (a) Transverse momentum 
distribution, (b) Spin alignment parameter a for the direct contribution, (c) Spin 
alignment parameter a with feed-down from \b decays taken into account. Dotted lines, the 
quark spin conservation hypothesis; dash-dotted lines, the full depolarization hypothesis. 
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